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Abstract 

The  pulsations  of  a gas-filled  bubble  in  a sound  field 
are  discussed  with  particular  attention  to  the  effect  of 
the  viscosity  of  the  external  liquid^  the  irreversible  con- 
duction of  heat  by  the  gas  within  the  bubbles  and  the  scat- 
tering of  energy  by  the  bubble.  The  expressions  obtained 
for  the  gas  pressure  within  the  bubble  and  the  time  rate 
of  change  of  bubble  radius  are  used  to  determine  the  rate 
at  which  gas  diffuses  in  and  out  of  a bubble  of  mean  radius 
R.  A threshold  for  bubble  growth.  Called  the  threshold  for 
gaseous^type  cavitation  is  defined  by  the  condition  that 
the  net  flow  of  gas  across  the  surface  of  the  bubble  will 
be  zero=  This  average  rate  of  gaseous  diffusion  can  also 
be  used  to  determine  in  a stepwise  fashion  the  mean  change 
of  bubble  radius  with  time=  Transient  effects  are  also 
discussed- 
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specific  heat  at  constant  pressure 
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Rotation  (continued) 

■ soles  of  gas 
= coshjrf  - eos*$ 

* molecular  weight 
« pressure 

■ pressure  tern  associated  with  viscosity 

■ gas  pressure  within  bubble 

■ variational  gas  pressure  within  bubble 

= scattered  pressure  wave 

= complex  amplitude  of  p’ 

s 

■ average  pressure  within  bubble  " p0  + 

■ hydrostatic  pressure 
= Legendre  polynomial 

a radial  particle  velocity  of  scattered  wave 
« heat  content 
*»  radial  coordinate 

■ radius  of  bubble 

« initial  radius  of  bubble 
= universal  gas  constant 
= time 

= initial  temperature  of  bubble  and  liquid 
“ interval  energy 
= volume 

* particle  velocity  vector 
“ work 
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Notation  (continued) 

Z ■ coaplex  iapedance  presented  by  bubble  to 

incident  sound  wave  = B + iV 

o 

Y * ratio  of  specific  hea'uJ  for  gas  ■ t* 

Gv 

5 ' /"f 

Y'l  X 

2/55“ 

0 - temperature 

0 * space  average  of  temperature 

K » wavelength 

p,  o viscosity 

p *■  density 

(T  - surface  tension 

/ <■  ySaHTE 

3 * scalar  potential 

to  * angular  frequency 

^ o resonant  angular  frequency  of  bubble 

In  general  the  subscripts  "1"  and  "2"  refer  respectively  to 
the  regions  internal  and  external  to  the  bubble. 
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I 

INTRODTTCTIOH 

This  report  presents  a theoretical  discussion  of  the  manner 
in  which  gas-filled  bubbles  in  liquids  react  to  the  incidence  of 
a sound  wave  and  then  grow  as  a result  of  a net  inward  diffusion 
of  gas.  This  process,  referred  to  as  rectified  diffusion,  re- 
sults when  more  gas  diffuses  into  the  bubble  during  its  half- 
cycle of  expansion  than  diffuses  out  of  the  bubble  during  its 
half-cycle  of  contraction.  The  theoretical  formulation  of  this 
problem,  namely  the  diffusion  of  gas  across  a moving  boundary, 
results  in  a nonlinear  equation.  As  a first  approximation,  the 
solution  to  the  diffusion  equation  is  obtained  for  a bubble  of 
constant  radius.  In  determining  the  rate  of  flow  of  gas  across 
the  surface  of  the  bubble  the  bubble  radius  (and  area)  is  then 
allowed  to  vary  in  accordance  with  the  pressure  variations  with- 
in the  bubble.  This  method  had  been  first  suggested  by  Harvey 
[1-6]*  who  considered  the  steady-state  diffusion  of  gas  into  a 
bubble  of  fixed  radius.  Harvey  assumed  that  the  liquid  surround- 
ing the  bubble  is  nonviscous  and  the  temperature  within  the  bubble 
constant.  Thir  work  was  later  extended  by  Blake  [7-10]  who  con- 
sidered the  periodic  steady-state  diffusion  of  gas  into  a bubble  of 
fixed  radius.  31ake  assumed  that  the  temperature  within  the  bubble 

♦References  to  articles  included  in  the  Bibliography  are  made 
by  numbers  enclosed  in  brackets  [ 3. 
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is  constant  and  that  all  of  the  energy  within  a sound  wave  inci- 
dent upon  the  bubble  is  used  to  cause  pressure  variations  within 
the  bubble,  Thus,  previous  work  with  regard  to  the  concept  of 
rectified  diffusion  had  neglected  the  effect  of  the  extinction 
of  energy  by  the  bubble  as  a result  of  heat  conduction  within 
it,  scattering  by  it,  and  the  resistive  forces  due  to  the  vis- 
cosity of  the  lj.quido  These  effects  are  discussed  in  this 
reports  In  addition,  the  transient  solution  of  the  diffusion 
process,  a factor  which  lias  not  been  examined  by  other  research- 
ers,  is  discussed,  In  order  to  obtain  an  expression  that 
describes  the  manner  by  which  gas  diffuses  into  the  bubble, 
relations  will  first  be  determined  for  the  variations  of  tem- 
perature within  the  bubble,  and  the  pressure  field  external 
and  internal  to  the  bubble « 


II 

TEMPERATURE  variations  within  the  bubble 


A,  The  Heat  Conduction  Equation 

Let  us  assume  that  a gas-filled  bubble  is  subject  to  a uni- 
form pressure  field  such  that  the  changes  in  pressure  at  internal 
points  are  due  to  the  fluctuations  of  its  boundary.  This  assump- 
tion holds  true  as  long  as  the  bubble's  diameter  is  much  less 
than  the  wavelength  of  the  incident  sound  wave.  During  the 
oscillations  of  the  bubble,  gas  will  diffuse  inward  and  outward, 
but  the  amount  of  gas  in  the  bubble  will  be  essentially  constant 
over  a single  cycle  at  practically  all  frequencies  of  interest 
in  this  report.  Within  the  bubble,  the  first  If.w  of  thermodynamics 
is  applicable  to  a small  element  of  volume,  v,  subject  to  a pres- 
sure, p:  Thus, 


dt  JFE 


dt 


* P 


(2-1) 


Schneider  [11]  has  shown  that  the  processes  of  evaporation  and 
condensation  are  sufficiently  rapid  in  comparison  with  the 
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pulsation  rate  of  the  bubble  so  that  the  pressure  of  the  vapor 
withtn  the  bubble  will  be  constant  and  determined  solely  by  the 
temperature  of  the  iiquido  In  general  the  pressure  of  the  gas 
within  the  bubble  will  far  excsed  the  pressure  of  the  vapor 0 
According  to  the  perfect  gas  law, 

pg  dv  + vdpg  • mg  R*  dOL,  (2-2) 


and  assuming  that  the  only  heat  conduction  process  takes  place 
In  the  gas 
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Then  In  spherical  coordinates  for  a spherically  symmetric  system. 
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where 


rO, 


The  heat  conduction  equation  external  to  the  bubble  is 


3^9  3u2 

ipr  * az  sr  * 0 


(2=7) 


* The  single  subscript,  1,  will  be  used  to  denote  the  region 
within  the  bubble?  the  single  subscript,  2,  will  be  used  to  de= 
note  the  region  external  to  the  bubble = 
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The  gas  pressure  within  the  bubble  will  be  aade  up  of  an  average 
pressure  tern  and  a variational  part  p^0  We  shall  see  in 
Chap»  III  that  we  can  write 

H • p’  •1“t  , 

and  the  coefficient  P|  nay  be  coaplexo 
We  shall  also  show  that 

?1  ■ Po  + ■ 


The  steady-state  solution  of  the  two  equations,  (2-6)  and 
(2-7) , can  be  obtained  by  writing  that 


Then 


u ■ v(r)es 


3 - *"Vl  - iieP^P  » 0 
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As  solutions  of  these  equations,  let  us  try 


■ A^inh  Q1r  + B1  cosh  q^r  - c (2-10) 
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v„  * A*  e + B-  e 
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Then,  since  v^ — >0  as  r — *0,  and  v2  is  finite  as  r — ,»od, 
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The  boundary  conditions  at  the  surface  of  the  bubble  based 
upon  the  continuity  of  teaperature  and  heat  flux  are 


e. 


Jr=B 
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lr=3 


Jr-B 


and  ^ ^(-^)  - K2  ^(-j2) 


Jr-R 


(2-15) 

r«R 


By  satisfying  these  boundary  conditions  we  find  that  the  teapera- 
ture  internal  to  the  bubble  is  given  by  the  expression 


[(l+OSgB+lJ  S Sinh(  1+1)83? 


(l+i)Cgj  S^coshd+OS^^^  sinh(l+i)81R3+(l«g^)sinh(l+i)51R 

P'  elwt  . (2=16) 
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This  solution  nay  be  written  as 
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P‘  eiut  , (2=17) 


where 


u0  ■ sinb5^r  ^osS^r  sinh  -fj  cos  + cosher  sinS^r  cosh  ^ sin  ^ 
vQ  * coshS^r  sin^r  sinh  ^ cos  ^ = sinh51r  cosS-^r  cosh  | sin  4 
U1 

U1  w" 
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The  ratio,  g » ^ , for  various 

licpiids  is  as  follows 

«!  - • i 

0, 055x10“-^ oal/seo  ca°C  for  air)g 
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Liquid 

(See  Saith  [12]) 

• 

i 1 

benzene 

0 ,14 

kerosene 

0,1? 

castor  oil 

0,13 

olive  oil 

0,14 

water 

0.04 

carbon  tetrachloride 

0,20 

acetone 

0,13 

When  r = R , u„  ■ a 1 , and  v * 0, 

’O  ’ 0 

Then  if  g « 1,  the 

tea- 

perature  at  the  surface  of  a bubble  (r  ■ R)  is  - TQ  0 = 
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A plot  of  u£«l  and  v£  as  functions  of  4 when  r » H is  presented 
in  Fig  = 1 and  Table  2 for  a given  value  of  g and  a'o 


Table  Z 


i 

4-1 

T1 

0,1 

- Q.OTxlO"3 

-0=l6zlC~3 

1 = 0 

- 2=  00xlC“3 

-2=13xl0“3 

*»10=40xl0“3 

-2=  55x10  ~3 

10o0 

-11=  80x10  ~3 

~l=37xl0“3 

IjOO.O 

-13.00xl0°3 

-0=Uxl0“3 

(g  = 0=20,  a® 

“ Va  2/al 

“ 15) 

Figure  1 represents  the  extreme  case  of  carbon  tetrachloride 
wherein  the  maximum  change  in  temperature  at  the  surface  of  the 
bubble  would  be  expectedo  It  is  apparent  that  even  if  the  peak 
pressure  within  the  gas  is  as  great  as  10  atmospheres,  the  sur- 
face temperature  will  differ  from  the  liquid  temperature  by  10=6°C 
(b-j/a-L  - =0o8l3xl0"4  for  air)  = 

By  means  of  these  curves  and  the  values  of  u0  and  vQ  one  can 
determine  the  temperature  at  various  points  within  the  bubble e 
For  convenience  one  can  let  ^r  - n ^ where  n is  any  positive 
fraction,  0 < n < 1,  such  that  R/r  ■ l/n0 

The  average  value  of  © with  respect  to  r will  be 

© ■ -4?  I r2e  dr.  (2-18) 

* 

By  using  the  integral 

^3  j r ” sinh  sr  dr  = cosh  sR  - ”2^2  sinil  sR’ 
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we  obtain  the  expression 
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where 


al  “ 


ga‘(cosh  if  -cos  rf)  + ■ | a'  if2(co3h  & + cos  tS) 


+ sinh  iftifa'^  a’+l~g)+$l-sini4[irfa,(^  a'+l+g)+l3 

f [J( 2 a'+l-g)2+  | a,23(cosh  cosrf)  + g2rf(cosh  £+  cos  if), 

+ 2g(rfa'+l-g)sinh  if  +2g(l-g)sin  if  j - 1, 

$1  ■ ^ /[(|  a'  + 1)(|  a1  + 1-g)  | + I a'2] (cos  if-  cos  t&) 

-grf(^  a'  + l)(oosh  6 + cos  rf) 

2 

+[2g<!  a'+l)  - (^~  a1 2+  ifa * 4 lflsinh  if 

2 

+[2g(|  a'+l)  - (^-  a'2+  rfa'  + l)3sin  if 

i^2  a'+l°g)2  + ^ a,23(cosh  if=cos  if)+g2if( co3h  if  + cosif) 

+ 2g(ifa'+l-g)  sinh  |f  + 2g(l-g)sinif  « 

A plot  of  a1  and  3^  is  presented  in  Pig0  2 and  Table  3 9 on  the 
opposite  pageo 

It  is  apparent  from  Pig,  2 that  ^->r  if  < 3»0  the  pulsations  of 
a bubble  will  be  essentially  isother^al^  for  if  > 30«0  the  pulsa- 
tions of  a bubble  will  be  essentially  adiabatiCe  In  the  transi- 
tion region  losses  will  occur  as  a result  of  the  irreversible 
conduction  of  heate  We  shall  see  later  that  the  resonant  radius 
of  a bubble  usually  lies  in  the  vicinity  of  this  transition 
region0 


FIG.  1.  THE  TEMPERATURE  AT  THE  SURFACE  OF'  A 


FIG.  2.  THE  SPACE  AVERAGE  OF  TEMPERATURE  WITHIN  A 


-9- 

Tafelfl  2 


O0  20 
15=0 


5=o 

10,0 

100,0 


xlO"3 

_ - G3^ 

xlO  - 

xlO“3 

x!0“3 

xlO' 

“3  xlO"3 

23  = 3 

- 1,0 

20,6 

0 

26,1 

0 

- 2,2 

- 33  = 7 

- 3 = 1 

- 34.8 

- 3 = 8 

- 35=3 

-391=8 

-354,0 

-395=5 

-353  = 0 

-397=3 

-352=5 

-701,0 

-239=4 

=702*9 

-243,1 

-703,4 

=238,0 

-970,1 

- 29.3 

-970.3 

- 29.1 

-970.4 

- 29  = 2 

Ihe  perfect  gas  law  for  the  bubble  is 

(Pi+  Pp  | ^R3  “ BgR*0, 


Now  let  R vary  with  Then 


.fJB}  dSL  . i] 

£**  dp*  J3(P± 


+ pi) 

Pg' 


But  by  differentiating  Bq,  (2-20), 
do?  ' V (al+  1S1> 


f ■ ^>-1  3^ 


The  quantities 
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p V 


Thus  Eqe  ( 20 23 ) becoaes 


4pi 


f * ‘V*  5(P4+  p'JT  - 


where 


Y - Sfi 

ST 


a * (y-Daj  + y 

3 ■ -Cy=1)3i 


(2-2*) 


This  equation  expresses  the  relation  between  the  change  of 
the  radius  of  the  bubble  and  the  pressure  variations  within  the 
gas-filled  bubbles  This  equation  will  be  used  later  to  study 
the  pressure  field  external  to  the  bubbles 


C„  Transient  Solution 

The  coaplete  solution  (transient  and  steady-state  terms)  to 
the  above  problem  can  be  obtained  by  means  of  Laplace  transforms 
(see  Carslaw  and  Jaeger  [133,  p»  288),  We  have  seen  that  the 
temperature  at  the  surface  of  a bubble  will,  in  most  cases,  not 
differ  greatly  from  the  temperature  of  the  surrounding  liquids 
Therefore,  the  importance  of  the  transient  term  can  be  obtained 
from  an  analysis  that  assumes  that  ths  fluid  reservoir  in  contact 
with  the  bubble  is  large  enough  and  has  a sufficiently  high  heat 
conductivity  such  that  its  tenperature  is  a constant  value  Tq0 
At  time  t » 0,  the  temperature  within  the  bubble  is  T0,  and 
the  pressure  within  the  bubble  is  PQ+  In  addition,  we  shall 

momentarily  assume  that  the  radius  of  the  bubble  is  constant  (see 
Appendix  I)0 

Let  us  define  the  Laplace  transform  of  ©,  namely, 

Q m f e“St  ©(r,t)dt  (2-2?) 

4 


TB2? 
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Equation  ( 2=5)  becoaes 

+ I -ii  - StsS,  + anO^  (t=o)-fc1  / e‘st  ££*  = A (2-26) 


■di  - a1s®1  4 a1®1ct“05-fc1  J ®'st  |£i 


dt  = 0 


d2(rSO 


CD  dp 

ll£r®l  + ,Air®i(  t^oj-b^r  ^ &~3*  dt  * 0 (2-27) 


The  initial  and  boundary  conditions  can  then  be  written  as 
0(t-o)  - T. 


®(r,s)  finite  as 
§(R,  s)  - Is 


But  since  Pg  “ Pg  “ pi 


dt  ■ f e”3^  tt®  = splc 


Thus  Eq„  (2-27)  becoaes 


-2(r§1)-a1sTl§1+  a-jTT^bjTsgJ  * 0 


(2-28) 


As  a solution  of  this  equation  let  us  try  the  substitution 


r®^  « Asinh  qr  + B cosh  qr  + Cr 


It  follows  that 


q - Va^s  , 


To  bl  - 

C ■ -A  - —*•  p'  . 

s a,  vg 


By  wirtue  of  the  boundary  conditions  we  obtain 
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and 


Thus, 


3 » 0 
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-II  fk  - 
al  3 sinh  Va^s  R 


5 . h.  f, 

a= 


rs  sinh  v'aTs  r 

-1 . i 


1 ° it  slnhva^s  B 


♦ V 


(2-29) 


The  Inverse  transfora  of  this  equation  is  obtained  by  aeons  of 
the  inversion  theorea 


f(t) 


. jl  r 

2ni  j 
Y*i  OD 


Y+i  ao  t 
" /*  fU) 


(2-30) 


The  tern  p^  is 


P' 


= s - ito 


By  choosing  the  proper  contour  in  the  conplex  plane  and  by  apply- 
ing the  theory  of  residues,  it  follows  that 
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b^VV  - Fi 0 


iwt 


+ p« 


”R  sinh  yiuaT  r 
_X-  -i 

_r  sinh  /iioa^  R 

jaSnSt 

2nn( -l)ne  aR^ 


(2-31) 


sin  nrr  | 


aR2(fl^a5  + Iw) 
aRz 


The  average  value  of  9 with  respect  to  r is 

K 

7 

9 - 


r29  dr, 


and  the  integrals 
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2 i r2  ^ sinh  qr  dr 


4?  cosh  qR 


2 2 sinil  qR 
Q « 


4?  f r2  ^ sin  nrr  § dr  ■ “(-1) 
R3  4 * R nrr 


n+1 


(2-32) 


a,B‘' 


In  the  majority  of  cases  the  tine  constant  -rg-  will  be  at 
most  of  the  order  of  magnitude  of  10 "A  seconds,  and  the  transient 
tern  can  be  neglected. 


Since 


■ ■ 3 ■ coth  vTjjaT  R = jj— - [sLnh/SjaT  R - sinysba,  R 

yiwiT?  R 1 /aua,  R 1 1 

1 1 (2-33) 

- i(sinhV2oa.,  R + sinV&a^  R)3, 

we  can  write  that 


where 


(a{+  ifipPj  eiijt 


And  as  we  would  expect  from  Eq.  (2-20),  the  terms 


(2-3*) 


1 


as  g- 


0, 
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III 

THE  PRESSURE  PIBLD  WITHIN  THE  BUBBLE 

Ao 

The  acoustical  field  external  and  internal  to  a spherical 
gas  bubble  can  be  determined  by  an  analysis  similar  to  one  devel- 
oped by  Epstein  [14]  3 The  small  signal  wave  equation  for  a vis- 
cous liquid  is 

p 0 +&xV*  §f  - fjVcV’  |f}-po |V(V.r,  . 0 (3-1) 

where 

•l  ■ <$>*  • 


By  writing  in  terms  of  a scalar  and  vector  potential 


v --v« 

+ v*r, 

(3-2) 

and  assuming  periodic  time  dependence,  it  follows  that 

+ k2  & - o, 

v2^  + ic2!*  - o , 

(3  “3) 

where 

k2  » w2[c2  + | i 

, E2«  -iw  ^ e 

Consider  the  case  of  a purely  longitudinal  primary  wave, 

A sinwt(Ia  Ae^*5),  that  i3  incident  upon  a gas  bubble  in  a vis- 
cous liquid.  If  the  origin  of  a spherical  coordinate  system  is 
taken  at  the  center  of  the  bubble  with  the  polar  axis  in  the 
direction  of  the  incident  wave,  axial  symmetry  will  require  that 
Ar»  Aq-  Oc  Both  9 and  A^  can  be  expanded  into  the  following 
series g 

= eikrcose  a ^ in(2n+i)jnac2r)Pn(cos0) 

2 n=o 


■14= 
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-1 5° 


S -^~V(  2n+l)B  hn(ksr)Pn(c030) 


n=o 

J3Q. 


V “ 

r2  n=o 


a> 


^>1°(2n-)CnA;Vi^  ?n(cosai 


*x  - \ in(V)PDCoos0) 

n=o  ^ 


where  the  subscript  "i"  refers  to  the  incident  wave$  the  subscript 
Mr’'  refers  to  the  reflected  wave$  the  number  "2"  refers  to  the  ex- 
terior of  the  bubble  $ the-  ^uaber  Mln  refers  to  the  interior  of 
the  bubble  $ and 


Vx) 


2 "n-** 


J-j-iCx)  o spherical  Bessel  function 


Vx)  * fS  Hn+^(x)  » spherical  Hankel  function 
PR(y)  “ Legendre  polynomial. 

Let  k2R  ■ c2,  kgR  * dg  an3  ^lR  " ci 5 Tlien  for  the  llcluids  and 

bubble  sizes  with  which  we  shall  be  concerned,  c?  « 1,  c|«l, 

k9  ■»  and  k-i  » » 

2 CL.  1 CL1 

4 1 

The  boundary  conditions  are  such  that  the  components  of  v’ as 
well  as  the  stresses  must  be  continuous  at  the  surface  of  the 
bubble  with  due  allowance  for  capillary  forces.  The  expansions 
for  Jn(ae)  and  h^z)  are? 


_n 

4_(x)  ■ — ~ — s - 

n lo3«  o » (2n+l) 


1 - 


2(2n+3)  2»<.(2n+3)  (2n+5) 


J_n-1(x)  - (-l)n  l83o«o(2ntl): 


»n~l 


- 

1 + 


a £ 0, 

(3 


L 2(2n-l)  2.4(2n-l)(2n-3) 

n > 1 


0125 
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J [i“|r  + fr  - 


and 


Vx)  = Jn(x)  " i^1)n^-n-l(x) 


Fop  the  first  tern  of  the  series,  that  is,  n ■ 0,  one  can 
by  using  these  expansions  for  J (x)  and  h^x)  that 

3 IB 

12  + * „2  = 

* o G “j  *“  C o 

Cg  0^  1 2 

<d i 


where 


Y = 1 - , f = 12 

* 2 3tPj_  ’ 


and  it  has  been  assumed  that 

pg  « pr 

Elimination  of  BQ  between  these  two  equations  leads  to 


B cs 
°1 


M'1- 

d2  J 


Hence  the  first  term  of  the  series  for  5^  is 


P2CL- 


y : Tzpn  + i ^ 


But 


and 


YP< 


plel1“'  o^TJ 


show 


(3=6) 


(3-7) 


(3-8) 


(3-9) 


(3-10) 
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Thus 


P ’ ~ 

r a ^ 


ct-i£ 


1 

P<Y-^  4 1 4&ZU 

3 1 3 


A e 


Iwt 


(3-H) 


The  coefficient  B,  can  be  determined  by  taking  the  boundary 
condition  for  n=l,  subject  to  the  same  approximations  used 


in  obtaining  B 


One  can  easily  show  that 

»!,-  % 


n ± j 2 

7 T ^2 


'1  7 

'1 


(3-12) 


The  additional  term  due  to  B,  in  the  series  expansion 


for  5^  is 


1 

7 + di 


Denote  the  correction  to  p^  due  to  this  additional  term  as 
Then 

-1 


P" 

Fg 


tf 

g 


»i 


r=R 


3 c2 


*1 


+ i 


ja L 


«L. 


(3-13) 


Thus  if  kR«l,  the  pressure  distribution  over  the  surface  of 
the  bubble  will  be  uniform*  For  bubbles  whose  size  becomes 
comparable  to  a wavelength,  the  additional  terms  of  the  series 
cannot  be  neglected. 


B,  Approximate  Analysis 


Equation  (3-11)  for  p^  can  be  obtained  in  far  more  simple 
manner  [153  if  one  makes  the  additional  assumption  that  the  terms 
pertaining  to  viscosity  in  the  wave  equation  are  negligible.  Then 
the  scalar  potential  for  spherically  symmetric  pulsations  in  the 
liquid  is  simply 


a - s0 


ei(wt  - kr) 

*r 


o 


(3°14) 


TM2  5 
Thus 


where 


p at 


*s  = 


Ps  ■ pressure  in  radially  scattered  sound  wave, 

q =.  particle  velocity  in  radially  scattered 
3 sound  wave. 

Up  to  this  point  the  effect  of  the  viscosity  of  the  fluid 
has  been  neglected.  As  a result  of  viscosity  an  additional 
pressure  tern,  p^,  will  appear.  The  approximate  value  of  pf  can 
be  obtained  in  the  following  manner0  The  pressure  tensor  in  a 
liquid  of  viscosity, u,  is,  for  the  case  of  radial  symmetry, 


Pf*^  a|Cr2qg)-  2u4f 


(3-16) 


The  pulsation  of  the  bubble  will  be  in  accordance  with  the  thermo- 
dynamic equation  (2-24), 


az  = a 


fi  zZ. 

Y 


pi  ’ 


where  it  has  been  assumed  that 


H <K  pi  • 

At  the  surface  of  the  bubble  the  continuity  of  pressure  and 
velocity  requires  that 


Pi+  pg 


Aeiajt  + p 


s]  + pf)  + ■¥" 

'r=B  / r=H 


p'  = -^r  (a  eiayt>+  p + p. 
dt  ^g  dt  \ vi 


IS  •££  (1-17) 
_2  dt 


and 
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(3-18) 


The  average  value,  Pi , with  respect  to  time,  is  simply  (PQ+  ^)  ? 
where  R is  a mean  value,  Thus  we  can  obtain  the  expression  for 
P’,  namely, 

3v  Pj/o-10  * t 
— * — A 


0 1 £$ 

- g 


■s-  p V 
J ? * 

a-ifT 


tu  ,^T»2 


STS  + 14uu  ♦ 1 


1+k  S 

which  is  the  same  as  Eq.  (3-11)  for  kR«  1. 


TT3 

fur- 5 

c(l+k2R2) 


C.  Impedance  of  Bubble  and  Resonance 


In  a similar  manner  we  can  determine  the  impedance  presented 
by  the  bubble  tc  a sound  wave,  namely, 

.iwt 


where 


Z - Mr 

(Vr=R 

- U + 17 


Ay-pfcfa ! 

R ^1+k2B2 


(3-20) 


7 « ttfiS 


2HiL_ 

wR(a2+02) 

1 


l-Hc^2  ^R(a2t_gS) 

3rPia-|i(a2+g2) 


The  termsin  U result,  respectively,  from  (1)  the  **ork  accomplished 
by  the  incident  pressure  in  overcoming  resistive  forces  due  to  the 
viscosity  of  the  liquid}  (2)  the  energy  scattered  by  the  bubble} 
(3)  the  energy  lost  as  a result  of  the  irreversible  conduction 
of  heat  within  the  bubble.  The  terms  in  V respectively  result 
from  (1)  the  inertial  reaction  upon  the  bubble  of  the  entrained 
fluid  that  moves  with  the  surface  of  the  bubble}  (2)  the  compres- 
sibility of  the  gas  within  the  bubble. 


A plot  of  the  terms  U and  V as  functions  of  ti  at  a frequency 
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of  60  kc/s  is  presented  in  Fig*  3 and  Table  4, below*  Two 
liquids  have  been  considered,  namely,  water  and  olive  oil0  The 
bubble  is  assumed  to  be  air-filledo 


Table  4 

Water  01iv_e_Ml 


J 

V 

S(  cm) 

u 

V 

tJ  V 

0=02 

1x10°^ 

400*0 

~8*52xl06 

336000  =4,31xl06 

0=1 

5*1x10“? 

790  o 5 

-4*59x10? 

66403  a 2 -2=94x10? 

1 

5dxio“4 

79,3 

-Ie86xl04 

6640=3  -Io71xl04 

5 

2<,5xl0“3 

64  * 2 

-2068xlO^ 

1374*7  -2*69xl03 

10 

5*  ixio~3 

50*4 

- 85*7 

704*6  - 216*9 

100 

5»ixio=2 

2530 

+I*89xl04 

2406  +1*  73*104 

The  steady-state  response  (or  impedance)  of  the  bubble  in- 
dicates to  some  extent  the  appearance  of  its  transient  pulsa- 
tions 0 The  impedance  consists  of  resistive,  inertial,  and 
stiffness  terms*  Let  us  assume  that  these  terms,  to  a first 
approximation,  are  constanto  For  small  values  of  H,  the  reac- 
tance will  be  negative  ==  that  is,  the  bubble  will  be  stiffness” 
controlled*  Under  these  conditions  one  can  show  [16]  that  it 
is  possible  for  an  underdamped  or  oscillatory  condition  to 
exist  in  "non-viscous"  liquids  wherein  a high=frequency  transient 
of  rather  large  magnitude  will  be  superimposed  on  the  steady-state 
pulsations  if  the  applied  signal  has  the  proper  phase = In  other 
words,  under  special  conditions  the  transient  pul ; .tions  of  a 
bubble  could  be  important  in  affecting  the  initial  diffusion  of 
gas  inward  and  outward.  We  shall  not,  however,  discuss  this 
problem  in  this  memorandum* 


Equation  (3-19)  can  now  be  written 


P L = 


“ IwRZ 


A e 


iwt 


P i JS 

' g 


lA) 


(3-21) 


TM25 


Resonant  pulsations  of  tbe  bubble  will  occur  when  V ■ 0,  or 


2/„2\o2  ’ 


l~(uo?/c*)R 


(3=22) 


where 


2 2 

007  " 0) 

1 0 


1 + £t-r  . ad &&) 

1 RPrtu  3r  ; 


2 S3  X 


* “o  B p2  p 


Similarly,  we  can  define  a resonant  bubble  size  R =>  Rf,  when 


oo  » coro  Thus , 


Z « U + i -^-§-s(l  - *4“) 


l+Jc2R2 


2 7 » 


(3  “23) 


and  the  displacement  of  the  surface  of  a bubble  is,  respectively, 


I 


r-R  /lu  3<V 
if- 1 — h- 


Similarly, 


-A  iftftJL. 

1 + k2R2 


3 


3Y<P0+  g^)/a-i$ 

iR~ 

,s  )r=>R 


for  R « R. 


for  R ■ R. 


for  R » R, 


for  R « R. 


for  R * R. 


3r<V  f) 
1R  3s)  _ 


for  R » R. 
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The  strain  set  up  at  the  surface  of  the  bubble  is  ~"r=R 

P 

This  strain  will  be  very  large*  The  corresponding 
pressure  will  not  be  very  largeo  However,  if  any  material  is 
near  or  in  contact  with  the  bubble,  it  may  be  mechanically 
damaged  as  a result  of  repeated  deformations  associated  with 
the  large  strains * Some  materials  will  effectively  change 
the  character  of  the  pulsations  of  the  bubble  such  that  the 
above  analysis  will  not  hold  true*  And,  finally,  the  effects 
of  numerous  bubbles  are  extremely  complicated  and  can  only  be 
estimated  on  the  basis  of  the  above  analysis* 


c 


IV 

THE  RECTIFIED  DIFFUSION  OF  OAS 

A*  Complete  Solution  to  the  Diffusion  Equation 

The  diffusion  of  gas  into  a bubble  will  be  in  accordance 
with  the  diffusion  equation 


For  a moving  boundary-value  problem  this  equation  is  nonlinear 
(see  Appendix  I)*  We  shall  assume  that  the  radius  of  the  bub- 
ble  is  momentarily  constant  so  that  an  exact  solution  to  this 
equation  can  be  obtained*  This  assumption  means  that  the  dif- 
fusion process  occurs  as  a result  of  variations  of  the  gas 
pressure  within  the  bubble.  These  variations  in  pressure  result 
from  the  continuity  of  pressure  across  the  surface  of  the  bubble* 
We  shall  later  malce  an  approximation  to  account  for  the  move- 
ments of  the  bubble 5s  surface* 

During  the  positive  half-cycle  gas  will  diffuse  into  the 
liquid  which  is  undersaturated  with  respect  to  the  bubble 5 during 
the  negative  half-cycle  gas  will  diffuse  into  the  bubble  which 
is  undersaturated  with  respect  to  the  liquid*  But  the  surface 
area  of  the  bubble  will  be  greater  during  the  negative  half- 
cycle such  that  there  will  be  a net  sonically  induced  diffusion 
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of  gas  into  the  bubble.  If  the  rate  of  this  net  influx  exceeds 
the  rate  at  which  gas  dissolves  as  a result  of  the  internal  ex= 
cess  pressure  due  to  surface  tension,  the  bubble  will  grow» 

This  process,  called  "rectified  diffusion,”  can  be  studied  by 
means  of  the  diffusion  equation^ 

In  spherical  coordinates  for  a spherically  symmetrie 
system  the  diffusion  equation  is 


A 

D 


(4-2) 


where 


The  boundary  and  initial  conditions  are 

c(r,t) — ► aQg0  as  r — ^ao 
c(R,t)  - aQpg 

e(r,o)  - a0g0  + § a0(pg(t=o)-g0), 


(4-3 ) 


or 


where 


u(r,t) vrg0  as  r — *•  cd 

u(R,t)  - Rpg  - H(Pj_+  a£A  sinCwt+D)  » W^t) 
u(r,o)  * rgQ  + H(pi"S0)  + R&jk  sinX1  3 Wg(r) 


X!  = h + i 

h = tan"1  y 

i ° tan'1  J 


j 


71125 


“21- 


The  oosplete  solution  of  the  diffusion  equation  in  the  infin- 
ite region  bounded  by  the  internal  sphere  of  radius  R is  readily 
obtainable  by  the  use  of  Fourier  integrals  and  Laplace  transforms* 
A complete  discussion  of  these  techniques  is  given  in  the  text- 
books of  Carslaw  and  Jaeger  [133  and  Churchill  [173*  The  solu= 
tion  is  (see  Carslaw  and  Jaeger  [133,  p*  209) 


<9 

- (r°A).B 

- lrtte2Rjf  " 

u - 1 

I W (X) 

e *Dt  - e 

4Dt 

2/r3)t 

K 2 

. 

- 

ah 


„ lr°Rl2  \Q=y2  dy 
4Dy2'  J 


(4"5) 


Let  * “ irik”  » 2 " i and  * (r+\-2R)f|.  Thens 


u 


Vt?” 


9, 


*tr-R)T) 

7 -i 

( r-S)Ti 


V.( 4 + r)  e-2  dz  - 
2V  R ' 


,r<40») 


W2(^  -r+2R)e“2  dz 


(r-R)2^ 
z2  * 


(4=6) 


e s dz 


where  the  integration  parameter  is  now  denoted  by  z*  Suppose  that 
the  boundary  conditions  are  those  of  Eq  = (4-4)*  Then  the  solution 
of  the  diffusion  equation  is 


u(r,t) 


rgQ+  R(Pi“g0)+  Ra^AslnX1 


erf 


r~R 

2^Dt 


+ Ra^A  e 


-Cr-R)  fM  , . . 

20  sinfwt  + X'-(r-R)^  jfe 


. ,n  r -Dtz* 

+ 2Ra*A  cosX 1 =*  / *§— - x-t 

1 n 6 w2+dV 


sinz(r-H)dz 


TK2? 
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- 2Ra{A  sinX'  ^ 


l2  T 

n J W2+D28A 


sinz(r-R)df = 


Thus , for  t > 0 , 


r»R 


*o(-s  *i  a) 


r=H 


* 3 


“ aoalA  V2D 


Ra4AsinX' 


P.,-  g * a, 'A  sin(o>t+X')  - — j; — 
1 0 1 VfrDt 

^ [sin(oit+X: ) +eos(wt+X  ' ) ] 

2 


+ 2a0a1A  cosX'^ 


9 „2  -Dtz£ 


,.0_ 


*4  + g4 


dz 


0 D2 


- 2a  ajA  sinX:  J f •M j 

ci  "J  u + s* 


The  term 


_2n  ~Dtz2 


dZ 


t > 0 


(4-7) 
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(*-9) 


has,  to  this  author's  Icnowledge,  not  previously  been  evaluated  in 
closed  form=  The  evaluation  will  therefore  be  given  in  consider- 
able detail  Quite  often,  integrals  of  the  form 


-az 


dz 


cannot  be  solved  by  Cauchy's  method  of  contour  integration  in 
the  complex  plane0  The  most  common  example  is  the  error  function 
of  infinity,  namely,  the  definite  integral 
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Other  techniques  in  the  real  plane  must  be  used  to  determine 
these  integrals,  The  integral  Rn  can  be  written: 


= jfi  vrar  f 

- S(a)t)  ) coscut- (-|  - C(wt)  )sinwt] , (4.-11) 


where  S(wt)  and  C(cot)  are  the  Fresnel  integrals  defined  by 

cot 

SCct)  - f -tig  d x 
J V 2nx 


co ' 


C(cot) 


cosx 

V 2ni 


dXc 


Thus  ? 
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C(^  - S(ut))sinwt  + (|  » C(wt))coswt] 

(4=»12) 


Similar  techniques  can  be  used  to  evaluate  many  integrals  of 
the  forms 


_2 

irf(*)e=az  da: 

The  associated  trigonometric  forms  of  these  integrals  may  also 
be  evaluated  by  means  of  integration  in  the  complex  plane « 


r 

A 


B.  Rectified  Diffusion 

Since  the  rate  at  which  gas  flows  into  a bubble  in  moles/sec 
is 

m - 4TtR2D(|£)  (4-13) 

r r-R 

we  can  now  determine  the  average  value  of  m over  one  cycle » How- 
ever, let  us  now  introduce  the  variation  R(t)« 

The  radius  R(t)  of  the  bubble  can  be  expanded  in  a Taylor's 
series  about  its  mean  or  initial  value,  RQo  This  series  is 


H * E + 

o ap> 


AD  + 

a 


o 


where 


or 


ps  ■ pi + H 


Ad  a p ’ , 
g 


Substitution  of  Eqso  (2-24)  and  (3=21)  leads  to  the  series 

/ 

A 


R = R0  1 + 


wR  v62+  V2 


sin(wt  + tan“^  + 
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HQC1  + a1A  3in(iot  + h)  r , 0 ) 


(4-14) 


Then 


m a =4  tt  Ba  D 
0 0 


pi“  g0+  ^i( cor^.  -sin-t  5(1.2H0/^  )A2] 


where  i.  * tan^/a, 

since  it  can  be  shown  graphically  that 


(4-15) 


S(a)t)sina3fc  ^ C(a)t)cosojt  C(ajt)eoa^t  s tract 
<%*  C(oct) sincat  «C(ut)cosajt  sin2cct- 
» C(ajt)sin^wt  5sS(ut)cos<jL)t 
^BS(a)t)ooswt  siaot 
«S  (cot)  co  soot  sin2ajt  « S(cot)sin\ot 
«0 
and 

S(cjot)sin2a/t  « 0c25 

C(o)t)sin2ojt  « 0„  25 

The  terns  in  Eq0  (4-16)  represent  the  effect  of  the  transient 
terms  in  the  solution  of  the  diffusion  equation  Even  if  the 
graphical  averaging  of  these  terns  is  taken  over  the  first  few 
cycles  of  the  applied  sound  signal,  the  average  of  these  terms  will 
be  approximately  rero0  Thus  the  transient  part  of  the  solution 
has  a negligible  effect  upon  the  average  rate  at  which  gas  moles 
diffuse  across  the  surface  of  the  bubble = 

In  order  to  have  (on  the  average)  no  moles  of  gas  entering 
or  leaving  the  bubble,  the  sound  amplitude,  A,  must  be  large 
enough  such  that  m = 0t  Let  us  denote  the  sound  amplitude  that 
satisfies  this  condition  as  the  threshold  for  the  growth  of 

gas-filled  bubbles  (also  referred  to  as  the  threshold  for  gaseous- 
type  cavitation) o Then, 


( 4“l6) 
(4=17) 
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A plot  of  A as  a function  of  j t is  presented  in  Fig = 4 and 
Table  5=  The  frequency  chosen  is  60  kc/s,  and  two  liquids, 
water  and  olive  oil,  have  been  considered,  The  bubble  is  assumed 
to  be  air- filled o 


Table  5 


i 

R(cm) 

Mel 

A ( Atmospheres) 

Olive  Oil 
A ( Atmospheres) 

0,02 

1*10=5 

14.8 

4 = 7 

0.1 

5=ixio"? 

1=9 

0.64 

1 

5*itio~4 

0o  14 

0 = 06 

5 

2.5xlO“3 

0c  02 

0.01 

10 

5oixio”3 

0,008 

0. 002 

100 

5=lxl0”2 

0,13 

0=05 

Co  Stepwise  Growth  of  Bubble 

For  a gas  bubble  nucleus  of  radius  R , much  smaller  than 
resonant  size,  a sound  amplitude  A greater  than  A ^ will  be  re- 
quired for  its  growth.  As  the  bubble  becomes  larger  the  sound 
amplitude  required  for  its  growth  (A  +)  becomes  less  and  less. 
If  the  initial  sound  amplitude,  A + , is  maintained  constant, 
the  rate  of  growth  of  the  bubble  will  be  more  and  more  rapid 
as  the  bubble  increases  in  size,  However,  because  of  the 
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additional  losses  that  occur  when  a bubble  approaches  its 
resonant  size,  the  rate  of  growth  will  be  retarded  for  the 
same  and  the  pulsations  of  the  bubble  will  become  very 

complex.  For  bubbles  greater  than  resonant  size,  the  rate 
of  growth  will  decrease  with  increasing  bubble  radius.  Hence, 
we  can  see  qualitatively  that  the  mean  radius  of  the  bubble 
should  follow  an  S-curve  with  time  for  a constant  value  of 
Naturally,  this  discussion  assumes  that  the  bubble 
remains  In  the  sound  field  at  all  times  and  is  not  subject  to 
any  forces.  In  a free  progressive-wave  system  this  assump- 
tion will  often  hold  true.  In  a standing-wave  system  Bjerknes 
forces  will  keep  a bubble  below  resonant  size  at  a pressure 
antinode  and  a bubble  above  resonant  size  at  a pressure  node. 
Because  of  these  complications  the  growth  curve  of  a gas- 
filled  bubble  can  be  best  obtained  by  means  of  numerical 
integration  of  the  nonlinear  equations  that  describe  the 
growth  of  the  bubble.  These  data  have  not  been  presented  in 
this  memorandom,  where  the  princip  points  of  interest  have 
been  (1)  the  minimum  threshold  for  growth  by  means  of  recti- 
fied diffusion,  and  (2)  the  characteristics  of  the  pulsations 
of  bubbles  of  various  sizes. 
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The  Nonlinearity  of  Parabolic  Equations 
Subject  to  Moving  Boundary  Conditions 


In  this  technical  memorandum  it  is  necessary  to  solve 
parabolic  equations  such  as  the  heat  conduction  equation  and 
the  diffusion  equation,,  For  these  equations  the  values  of 
the  temperature  and  gas  concentration  at  the  moving  surface 
of  the  bubble  can  be  taken  as  the  boundary  conditions,  A 
parabolic  equation  subject  to  conditions  at  a moving  boundary 
is  nonlinear  [18], 


Consider  the  one-dimensional  equation 


J L * - £1 

ax  F ax  " at  ♦ 


where 


or 


f(x,o)  = gCx)  , 


» 


R0<  S<  B!, 


x = R; 


f (R : , t)=  0 I 

1 J 

f(RCt)  ,t)=  h(t)  , x = R(t), 

Let  us  introduce  the  transformation 


(A-l) 


t B 0 


t > 0 
t > 0 


- R ! c (A~2^ 

Then 

1 S L £I_\  „ 3f  + JLmXj  - MU)  if-  0 < X ; < 1 , 
(R(t)-RO2  V 3XV  St  R(t)-R=  « ^ ’ t>0> 

(a-35 


where 
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f(x,o)  - g(x') 


II 

dx ' 


0 1 


xJ»l 


or 


f(l,  0)  - 0 j 
f(B(t) ,t)  - h(t) 


t » 0 


» t _ 


t > 0 


A.  Nfc. 


0. 


This  equation  is  nonlinear  unless  R(t)  is  constant*  There- 
fore, in  studying  the  pulsations  and  growth  of  gas  bubbles  we 
consider  a spherical  bubble  immersed  in  an  infinite  liquids 
This  bubble  is  in  the  path  of  a longitudinal  sound  wave* 

In  solving  equations  of  the  parabolic  type  we  find  it  conven- 
ient in  many  instances  to  assume  that  the  radius  of  the  bubble 
is  momentarily  fixedc  Later,  we  can  introduce  approximations 
for  the  variations  of  bubble  radius  with  tine,  namely,  R(t), 
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